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Abstract 

We exhibit some explicit co-adapted couplings for rt-dimensional Brown- 
ian motion and all its Levy stochastic areas. In the two-dimensional case 
we show how to derive exact asymptotics for the coupling time under 
various mixed coupling strategies, using Dufresne's formula for the dis- 
tribution of exponential functionals of Brownian motion. This yields 
quantitative asymptotics for the distributions of random times required 
for certain simultaneous couplings of stochastic area and Brownian mo- 
tion. The approach also applies to higher dimensions, but will then lead 
to upper and lower bounds rather than exact asymptotics. 
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Introduction 

It is a pleasure to present this paper as a homage to my DPhil supervisor 
John Kingman, in grateful acknowledgement of the formative period 
which I spent as his research student at Oxford, which launched me into 
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a deeply satisfying exploration of the world of mathematical research. It 
seems fitting in this paper to present an overview of a particular aspect of 
probabilistic coupling theory which has fascinated me for a considerable 
time; given that one can couple two copies of a random process, when 
can one in addition couple other associated functionals of the processes? 
How far can one go? 

Motivations for this question include: the sheer intellectual curios- 
ity of discovering just how far one can push the notion of probabilistic 
coupling; the consideration that coupling has established itself as an 
extremely powerful tool in probability theory and therefore that any 
increase in its scope is of potential significance; and the thought that 
the challenge of coupling in extreme circumstances may produce new 
paradigms in coupling to complement that of the classic reflection coup- 
ling. 

It has been known since the 1970s that in principle one can couple 
two random processes maximally; at first encounter this fact continues 
to delight and surprise researchers. I summarize this point in Section 
[T] and also describe the important class of co-adapted couplings. These 
satisfy more restrictive requirements than maximal couplings, are typic- 
ally less e fficien t, but are also typically much easier to construct. Since 
Lindvalll (1982)'s seminal preprint we have known how to couple Eu- 



clidean Brownian motion using a simple reflection argument in a way 
which (most unusually) is both maximal and co-adapted, and this has 
led to many significant developments and generalizations, some of which 
are briefly sketched in Section [5] This leads to Section [31 which develops 
the main content of the paper; what can we now say about the question, 
how to couple not just Brownian motion, say, but also associated path 
integrals? Of course we then need to vary our strategy, using not just 
reflection coupling but also so-called synchronous coupling (in which the 
two processes move in parallel), and even rotation co upling, which co r- 



relates different coordinates of the two processes. In iKendalll (|2007f ) I 
showed how to couple (co-adaptively) n-dimensional Brownian motion 
and all its stochastic areas; this work is reviewed in Section 13.31 using 
a rather more explicit coupling strategy and then new computations 
are introduced (in Section I3.4[) which establish explicit asymptotics for 
the coupling time for suitable coupling strategies in the two-dimensional 
case, and which can be used to derive naive bounds in higher dimensions. 
Section H concludes the paper with some indications of future research 
directions. 
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1 Different kinds of couplings 

Probabilistic couplings are used in many different ways: couplings (real- 
izations of two random processes on the same probability space) can be 
constructed so as to arrange any of the following: 



for the two processes to agree after some random time (which is to say, 
for the couplin g to be successful) . This follows the pioneering work of 



Doeblinl (|1938), which uses this idea to provide a coupling proof of 



convergence to equilibrium for finite-state ergodic Markov chains; 

for the two processes to be interrelated by some monotonicity property 
— a c ommon use of coupling in the study of interacting particle sys- 
tems (|Liggettl . [2005h : 

for one process to be linked to the other so as to provide some in- 
formative and fr uitful representation, as in the case of the coalescent 
(|KingmanL[l982h : 

for one of the processes to be an illuminatin g approximation to t he 
other; this appears in an unexpected way in iBarbour et al.l (|19920 's 
approach to Stein-Chen approximation. 



These considerations often overlap. Aiming for successful coupling has 
historical precedence and is in some sense thematic for coupling theory, 
and we will focus on this task here. 



1.1 Maximal couplings 

The first natural question is, how fast can coupling occur? There is 
a remarkable and satisfying answer, namely that one can in principle 
construct a coupling which is maximal in the sense that it maximizes 
the probabil it y of coupling b e fore time t for al l possible t: see lGriffeatb 
( 1974 Il978l ). IPitmanl |l976h . iGoldsteinl |l978l ). Briefly, maximal coup- 
lings convert the famous Aldous inequality (the probability of coupling 
is bounded above by a simple multiple of the total variation between 
distributions) into an equality. Constructions of maximal couplings are 
typically rather involved, and in general may be expected to involve de- 
manding potential-theoretic questions quite as challenging as any prob- 
lem which the coupling might be supposed to solve. Pitman's approach 
is perhaps the most direct, involving explicit construction of suitable 
time-reversed Markov chains. 
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1.2 Co-adapted coupling 

Maximal couplings being generally hard to construct, it is useful to con- 
sider couplings which are stricter in the sense of requiring the coupled 
processes both to be adapted to the same filtration. Terminology in 
the literature varies: Markovian, when the coupled processes are jointl y 



Markov, with prescribed m arginal kernels (jBurdzv and Kendall 120001: 



co-immersed (jEmervl . 120051 1 or co-adapted to emphasize the role of the 
filtration. The idea of a co-adapted coupling is simple enough, though 
its exact mathematical description is somewhat tedious: here we indicate 
the definition for Markov processes. 

Definition 1.1 Suppose one is given two continuous-time Markov pro- 
cesses X^ and X^ 2 \ with corresponding semigroup kernels defined for 
bounded measurable functions / by 

P®f(z) = E[f(xi%) | = z, X« for u<s]. 

A co-adapted coupling of X^- 1 ' and X^ is a pair of random processes 
X^ and X (2) defined on the same filtered probability space (f2, {3t : 
t > 0},P), both adapted to the common filtration {$ t ■ t > 0} (hence 
'co-adapted') and satisfying 

Pf/W = E[f(X$ t )\3.,X®=z] 

for i = 1, 2, for each bounded measurable function /, each z, each s, 
t > 0. 

Thus the individual stochastic dynamics of each JfW agree with those 
of the corresponding X^' even when the past behaviour of the oppos- 
ite process is also taken into account. (This is typically not the case 
for maximal couplings.) In particular, if the X^ are Brownian motions 
then their forward increments are independent of the past given by the 
filtration. Moreover if the processes are specified using stochastic dif- 
ferential equations driven by Brownian motion then general co-adapted 
couplings can be represented using stochastic calculus ( possibly at th e 
price of enriching the filtration) , as observed i n passi ng bv lEmervl (|2005l) , 



and as described more formally in iKendalll (|2009al Lemma 6). Briefly, 



any co-adapted coupling of vector-valued Brownian motions A and B_ 
can be represented by expressing A as a stochastic integral with respect 
to S and perhaps another independent Brownian motion C_: we use this 
later at Equation (|3.ip . 
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1.3 Coupling at different times 

There ar e many oth e r usef ul couplings falling outside this framework: for 
example iThorissonl (|1994l ) discusses the idea of a shift-coupling, which 
weakens the co upling requirem ent by permitting processes to couple at 
different times; iKendalll (|1994f ) uses co-adapted coupling of time-changed 
processes as part of an exploration of regularity for harmonic maps. 
However in this paper we will focus on co-adapted couplings. 



2 Reflection coupling 

The dominant example of coupling is reflec tion coup l ing fo r Euclidean 



Brownian motions A and B_, dating back to iLindvalll (|1982l )'s preprint: 
construct A from B_ by reflecting B_ in the line segment running from B_ 
to A. That this is a maximal coupling follows from an easy computation 
involving the reflection principle. It can be expressed as a co-adapted 
coupling; the Brownian increment for A is derived from that of B by a 
reflection in the line segment running from B_ to A. Many modifications 
of the reflection coupling have been derived to cover various situations; 
we provide a quick survey in the remainder of this section. 



2.1 Maximality and non-maximality 



The reflection coupling is u nusual in being both co-adapted and max- 
imal. iHsu and Sturml |2003) point out that reflection coupling fails to be 
maximal even for Euclidean Brownian motio n if the Brownian motion 
is stopped on exit from a prescribed domain. (Kuwad a and Sturmll2007 



discuss the manifold case; see also iKuwadal 120071 [2009J.) Perhaps the 
simplest example of an instance where no co-adapted couplin g can be 
maxi mal arises in the case of the Ornstein-Uhlenbeck process ([Connor , 



2007, PhD thesis, Theorem 3.15). Consider the problem of constructing 



successful co-adapted couplings between (i) an Ornstein-Uhlenbeck pro- 
cess begun at 0, and (ii) an Ornstein-Uhlenbeck process run in statistical 
equilibrium. A direct argument shows that no such co-adapted coupling 
can be maximal; however in this case reflection coupling is maximal 
amongst all co-adapted couplings. (The study of couplings which are 
maximal in the class of co-adapted couplings promises to be an inter- 
esting field: the case of random walk on the hypercube is treated by 
Connor and Jacka 20081 ) 
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2.2 Coupling for Diffusions 



A variant of reflecti on coupling for elliptic diffu sions with smo oth coeffi- 
cients is discussed in lLindvall and Rogers! (|1986[ ) and further in lChen and Li 

Reflection' here depends on interaction between the two diffu- 
sion matrices, and in general the two coupled diffusions do not play 
symmetrical roles. In the case of Brownian motion on a manifold one 
can use the mechanisms of stochastic development and stochastic par- 
allel transport to define co-adapted couplings in a more symmetrical 
manner. The behaviour of general co-adapted Brownian coupl ings on 
Riemannian manifolds is related to curvature. Kendalll |l986b[) shows 
that successful co-adapted coupling can never be almost-surely success- 
ful in the case of a simply-connected manifold with negative curvatures 
bounded away from zero. On the other hand a geometric variant of re- 
flection coupling known as mirror coupling will always be almo st-surely 
successful if the man i fold has non- negative Ricci cu r vatur es (jKendall 
1986al ICranstonlll992l lKendal]||l988h . Indeed iRenessel |20o4 shows how 



to generalize mirror coupling even to non-manifold contexts. 



3 Coupling more than one feature of the process 

The particular focus of this paper concerns ongoing work on the following 
question: is it possible co-adaptively to couple more than one feature of 
a random process at once? To be explicit, is it possible to couple not 
just the location but also some functional of the path? 

On the face of it, this presents an intimidating challenge: control of 
difference of path functionals by coupling is necessarily indirect and it 
is possible that all attempts to control the discrepancy between path 
functionals will inevitably jeopardize coupling of the process itself. 

Further thought shows that it is sensible to confine attention to cases 
where the process together with its functional form a hypoelliptic dif- 
fusion, since in such cases the Hormander regularity theorem guaran- 
tees existence of a density, and this in turn shows that gen eral coupling 
(not necessarily non-co-adapted) is possible in principle. ( Hairer 20021 



uses this approach to produce non-co-adapted couplings, using careful 
analytic estim ates and a regeneration argument which corresponds to 



Lindvaljl2002} s '7 coupling'. 

Furthermore it is then natural to restrict attention to diffusions with 
nilpotent group symmetries, where one may hope most easily to discover 



Distribution asymptotics for Levy stochastic area couplings 



7 



co-adapted couplings which will be susceptible to extensive generaliza- 
tion, parallelling the generalizations of the Euclidean reflection coupling 
which have been described briefly in Section [2J 



3.1 Kolmogorov diffusion 

Consider the so-called Kolmogorov diffusion: scalar Brownian motion B 
plus its time integral J B dt. This determines a simple nilpotent diffu- 
sion, and in fact it can be coupled co-adaptively by varying sequentially 
between reflection coupling and synchronous coupling fallowing the two 



Brown i an motions to move in para llel) as shown in iBen Arous et al 



( 19951) . Jansons and Metcalfe! ( 2007h describe some numerical investiga- 



tions concerned with optimizing an exponential moment of the coupling 
time. 

The idea underlying this coupling is rather simple. Suppose that we 
wish to couple (B,J Bdt) with (B,J Bdt). Set U = B - B and V = 
J B dt — J B dt. Co-adapted couplings include stochastic integral repres- 
entations such as dB = J dB, for co-adapted J E { — 1,1}; J = 1 yields 
synchronous coupling and J = — 1 yields reflected coupling. Suppose 
Uo 7^ and V = (which can always be achieved by starting with a 
little reflected or synchronous coupling unless U — V — from the start, 
in which case nothing needs to be done). We can cause (U, V) to wind re- 
peatedly around (0, 0) in ever smaller loops as follows: first use reflection 
coupling till U hits —Uq/2, then synchronous coupling till V hits 0, then 
repeat the procedure but substituting —Uq/2 for Uq. A Boref-Cantclli 
argument combined with Brownian scaling shows that (U, V) then hits 
(0 , 0) in finite time. 

Kendall and Price] (2004) present a cleaned- up version of this argu- 



ment (together with an extension to deal in addition with J J B ds dt). 

Curiously, this apparently artificial example can actually be applied to 
the study o f the t ail cr-algebra of a certain relativistic diffusion discussed 
bv lBailleull (l2008h . 



Remarkably it is possible to do very much better by using a com- 
pletely different and implicit method: one can couple not just the time 
integral, but also any fini te number of additional iterated time integrals 



(jKendall and Pried . 120041 ) , by concatenating reflection and synchronous 
couplings using the celebrated Morse-Thue binary sequence 0110100110010110 . 
Scaled iterations of state-dependent perturbations of the resulting con- 
catenation of couplings can be used to deliver coupling in a finite time; 
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the perturbed Morse-Thue sequences encode indirect controls of higher- 
order iterated integrals. 



3.2 Levy stochastic areas 

Moving from scalar to planar Brownian motion, the natural question is 
now whether one can co-adaptively couple the nilpotent diffusion formed 
by Brownian motion (B%, B2) and the Levy stochastic area J {B\ AB2 — 
B2&B1). This corresponds to c oupling a hypoell iptic Brownian motion 
on the Heisenberg group, and iBen Arous et al.l determine an explicit 
successful coupling based on extensive explorations using computer al- 
gebra. 

Again one can do better ( Kendall! . 2007 ). Not only can one construct 
a simplified coupling for the 2-dimensional case based only on reflec- 
tion and synchronous couplings (switching from reflection to synchron- 
ous coupling whenever a geometric difference of the stochastic areas ex- 
ceeds a fixed multiple of the squared distance between the two coupled 
Brownian motions), but also one can successfully couple n-dimensional 
Brownian motion plus a (T) basis of the various stochastic areas. In the 
remainder of this paper we will indicate the method used, which moves 
beyond the use of reflection and synchronous couplings to involve ro- 
tation couplings as well (in which coordinates of one of the Brownian 
motions can be correlated to quite different coordinates of the other). 



3.3 Explicit strategies for coupling Levy stochastic area 



Here we describe a variant coupling strategy for the n- dimensio n al cas e 
which is more e xplicit t h an th e strategy proposed in iKendalll (|2007T) . 
As described in iKendaill (|2007l Lemma 6), suppose that A and B_ are 
co-adaptively co upled n-dimensional Brownian motions. Arguing as in 
Kendalll (|2009al ). and enriching the filtration if necessary, we may rep- 
resent any such coupling in terms of a further n-dimensional Brownian 
motion C_, independent of B_; 



AA 



J T dB 



J T dC. 



(3.1) 



where J, J are predictable (n x n)-matrix-valued processes satisfying 
the constraint 



I 1 + I I = I 
and where I represents the (n x n) identity matrix. 



(3.2) 
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Note that the condition p. 21) is equivalent to the following set of 
symmetric matrix inequalities for the co-adapted process J (interpreted 
in a spectral sense): 

£ 1 < I- (3-3) 

Thus a legitimate coupling control ,]_ must take values in a compact 
convex set ofnxn matrices defined by (|3.3p . 

The matrix process J_ measures the correlation (dB_dA v )/dt 
between the Brownian differentials dA and dB_; for convenience let 
S_= h(J_ + jj) and A = |(J — J T ) refer to the symm etric and skew - 



symmetric parts of J. The coupling problem solved in IKendalll (120071) 



is to choose an adapted J = S_ + A which brings A and B into agree- 
ment at a coupling time T coup ij ng which is simultaneously a coupling 
time for all the Q) corresponding pairs of stochastic area integrals 
/ (Ai dAj - Aj dAi) and / (B t dBj - Bj dBi). 

To measure progress towards this simultaneous coupling, set X = A — 
B_ and define 21 to be a skew-symmetric matrix of geometric differences 
between stochastic areas with (i, j) th entry 

%ij = J (Ai d. l , - Aj dA t ) - J (Bi dBj - Bj dB t ) + A, 11, - A, 11,. 

(3.4) 

The nonlinear correction term AiBj — AjBi is important because it 
converts 21^ into a geometrically natural quantity, invariant under shifts 
of coordinate system, and also because it supplies a very useful contribu- 
tion to the drift in the subsequent Ito analysis. Of course 21 and X_ both 
vanish at a given time t if and only if at that time both A = B_ (so in 
particular the correction term vanishes) and also all the corresponding 
stochastic areas agree. 

Some detailed Ito calculus (originally carried out in an impleme nt a- 



2001, but 



tion of the Ito calculus procedures Itovsn3 in Axiom, IKendalll 
now checked comprehensively by hand) can now be used to derive the 
following system of stochastic differential equations for the squared dis- 
tance V 2 = \\2£}\ 2 and the 'squared areal difference' U 2 = trace(2J. T 2l) = 

(d(V 2 )) 2 = 8 E T (l-g)EV 2 dt, 

Drift d(V 2 ) = 2 trace (l - dt, 

d(V 2 ) x d(U 2 ) = -lQvJ^jA z vUV 2 dt 1 

( d(U 2 )) 2 = 32 vJZl + ZlK U 2 V 2 dt, 
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Drift d(U 2 ) = 4trace(j[ T ^) U dt + 

+ 4 (trace (| + g) - vj (| + g) v) V 2 dt . 

(3.5) 

Here the vector v_ and the matrix Z_ encode relevant underlying geometry: 
respectively Vv_ = X_ and UZ_ = 21. Note that y_ is a unit vector and Z_ 
has unit Hilbert-Schmidt norm: trace Zj Z_ = !• 

The strategy is to consider V 2 and £/ 2 on a log-scale: further stochastic 
calculus together with suitable choice of coupling control J then permits 
comparison to two Brownian motions with constant negative drift in a 
new time-scale, and a stochastic-calculus argument shows that K = 
^log(V 2 ) and H — |log([/ 2 ) reach — oo at finite time in the original 
time-scale measured by t. In fact further stochastic calculus, based on 
the martingale differential identity 

dlogZ = ^-\{^r " 
8 Z 2 \ Z 

shows that 

{dKf = \vJil-gvdT, 

Drift (dK) = \ (trace (| - g) - 2 vj (l - g) v) dr 

dK x dH = - E T g T AE w ' 

(dH) 2 = 2 v T g T {l + ggv d? 

Drift (dH) = \ trace (gJAjw 



W 1 ' 



+ | (trace (| + D - i/ T (1 + g) v - 4 i/ T ^ T (l + £v\ 



dr 

U -2 



(3.6) 

Here dr = 4w| determines the new time-scale r, and W = U/V 2 — 
exp (iJ — 2K). It is clear from the system (|3.6I) that the contribution of 
| trace (^ T A) w ( d 

eriving ultimately from the areal difference correc- 
tion term in (13.41) 1 is potentially a flexible and effective component of 
the control if 1/W 2 is small. On the other hand if this is to be useful 
then I — S_ must be correspondingly reduced so as to ensure that H and 
K are subject to dynamics on comparable time-scales. 

Kendalll ( 2007 ) considers the effect of a control J which is an affine 



mixture of reflection I — 2v_vJ and rotation exp (— 9 2±) . Second-order 
Taylor series expansion of the matrix exponential is used to overcome 
analytical complexities at the price of some mild asymptotic analysis. 
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Here we indicate an alternative route, replacing the second-order trun- 
cated expansion J' =\ — 0Z_ — \d 2 Zj Z_ (which itself fails to satisfy Q3.3P 
and is not therefore a valid coupling control) by J" = I — Z_ — 9 2 Zj Z_, 
which does satisfy (|3.3I) when 9 2 < 2 (using the fact that non-zero ei- 
genvalues of Zj Z_ have multiplicity 2, so that Trace(Z T Z) = 1 implies 
that g < gjg < A I). Thus we can consider 

1 = I- 2PKK T - (1 -p)9g- (1 -p)0 2 £ T £, (3.7) 

which is a valid coupling control satisfying (|3.3[) when < p < 1 and 
6 2 < 2. 

This leads to the following stochastic differential system, where terms 
which will eventually be negligible have been separated out: 

(dK) 2 = (p+i||i^||V) dr - l\\gy\\ 2 e 2 At, 

Drift(dif) = -\ (p - ± (l - 2\\gv\\ 2 ) 6 2 ) df 

- I(l-2||^|| 2 )0 2 df, 

dKxdH = \\gufe -f - p\\gvfO # , 

(dif) 2 = 4||^ £ || 2 ^ - 2(l-p)||^ T ^|| 2 2 ^, 

Drift(dtf) = -i0f + (n-l-4||^|| 2 ) ^ 

+ hp w - V (i - U4 2 - 4 iii T i^ii 2 ) o 2 ■ 

(3.8) 

In order to ensure comparable dynamics for K and H, set p = a 2 /W 2 
and = P/W (valid when > max{-^/3, a}): writing dr = dr/W 2 = 
A(V/U) 2 dt leads to 

(dK) 2 
Drift (dK") 

d-K" x dH 

(dH) 2 
Drift (dff) 



(a 2 + ^||ii£|| 2 /3 2 - ^li^ll 2 ^) dr, 
-Ml dr , 

lli^| 2 (l-^)/3 dr, 

4 (|li^|| 2 -|(l-^)|li T id| 2 ^) dr, 
-ua dr, (3.9) 
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where jjbx and \i2 are given by 

Ml = ia 2 -i(l-2|li£|| 2 )/3 2 + 3 ^(l-2||^|| 2 )a 2 /? 2 , 
fi 2 = \P - n + 1 + A\\Zlv\\ 2 

- 2m - (1 - £) (1 - iii^ii 2 - 4 iii T i-ir) p 2 ) ■ 

(3.10) 

In order to fulfil the underlying strategy, /ii and M2 should be chosen 
to accomplish the following: 

1. W = exp(H — 2K) must remain large; this follows by a strong-law- 
of-large-numbers argument if fj,\, fi 2 are chosen so that 2/ii — fi2 = 
Drift d(H — 2K) / dr is positive and bounded away from zero and yet 
(AH) 2 / dr and (dK) 2 / dr are bounded. 

2. Both of -/ii = Drift d(K)/dr and -fj, 2 = Drift d(iJ)/dr must re- 
main negative and bounded away from zero, so that K and H both 
tend to —00 as r — > 00. 

3. If coupling is to happen in finite time on the i-time-scale then 

1 coupling - 4 



exp (2H - 2K) dr < 00 . 

This follows almost surely if fii, \ii are chosen so that 2/ii — 2\xi = 
Drift d(2iJ — 2K) / dr is negative and bounded away from zero. 

Now the inverse function theorem can be applied to show that for any 
constant L > there is a constant V > such that if W 2 > V then 
(|3.10[) can be solved for any prescribed < fix, fJ-2 < L using a 2 , (3 2 < L' 
(incidentally thus bounding (d_ff) 2 /dr and (dK) 2 / dr). Moreover by 
choosing L' large enough it then follows that p — a 2 /W 2 < 1 and 
f3 2 /W 2 < 2, so that the desired /zi, 112 can be attained using a valid 
coupling control. 

A comparison with Brownian motion of constant drift now shows that 
if initially W > 2L' then there is a positive chance that W > L' for all 
r, and thus that coupling happens at r-time infinity, and actual time 
^coupling < 00 ■ Should W drop below L' , then one can switch to the 
pure reflection control (p — 1) and run this control until W = U/V 2 
rises again to level 2L' . This is almost sure to happen eventually, since 
otherwise V —> and thus U = WV 2 —> 0, which can be shown to have 
probability under this control. 
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3.4 Estimates for coupling time distribution 



In the planar case n = 2 the stochastic differential system (J3T6J) un- 
der mixed rotation-reflection controls can be simplified substant ially. 
In this s ection we go beyond the work of iBen Arous et al and 



Kendalll (|2007f) by using this simplification to identify limiting distri- 
butions for suitable scalings of the coupling time T coup ii ng . The simpli- 
fication arises because the non-zero eigenvalues of the skew-symmetric 
matrix Z_ have even multiplicity; consequently in the two-dimensional 

case we may deduce from trace(Z T Z) = 1 and = 1 that ||i?£|| 2 = \ 
and \\Zj Z_v_\\ 2 ~ j. Moreover an Euler formula follows from Z_Z_ = — 5I: 

exp (-V29g\ = cos0|-\/2 sin Z . 
Accordingly in the planar case the mixed coupling control 

£ = p(\-2vvj) +qexp(-V26£) 

= I- 2-pvvJ - q (1 - cos 0)1- V2 q sin Z_ (3.11) 

(for 0<p=l — g < 1, unrestricted 0) renders the system p.6p as 

( dK) 2 = (p+ \q{l - cos0) dr, 
Drift (dK) = -ip dr, 
dKxdH = -i=gsin0^, 

(dH) 2 = (2 - g(l- cos 0)) 
Drift(dif) = - 4= gsin0f -(1- | g (l-cos0)) & . (3.12) 



If we set = (so that the coupling control is a mixture of reflection 
coupling and synchronous couplings) then the result can be made to 
yield a successful coupling: choosing p = min{l, a 2 /W 2 } the stochastic 
differential system becomes 



(dK) 2 = mm{W 2 ,a 2 } dr , 

2 



Drift (d2iQ = -±min{W 2 ,a 2 } dr , 

dK x dH = , 

(dH) 2 = 2 dr, 

Drift(diT) = -dr, (3.13) 

where dr = Adt/V 2 as before. Accordingly, for fixed a 2 , once W 2 > a 2 
then H and K behave as uncorrelated Brownian motions with constant 
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negative drifts in the r time-scale; moreover if a 2 > 1 then 

Drift d(log W)/ dr = a 2 - I 

is strictly positive and so W — ¥ oo almost surely. In order to argue as 
before we must finally show almost-sure finiteness of 



T, 



coupling 



1 r°° 

- J exp (2H - 2K) dr . 



Now if we scale using the ratio (Uq/Vo) 2 at time zero then we can deduce 
the following convergence-in-distribution result as Wo = Uq/V 2 — > oo: 



U 



coupling 



i r°° 

- / exp (2H -2K - {2H - 2K )) dr 
4 Jo 

exp (W2 + a 2 B — (2 



a 2 )r 



dr (3.14) 



where B is a standard real Brownian motion begun at 0. This integral 
is finite when a 2 < 2, so finally we deduce that coupling occurs in finite 
time for this simple mixture of reflection and synchronous coupling if we 
choose 1 < a 2 < 2. 

However we can now say much more, since the stochastic integral in 
(|3.14[) is one of the celeb rated and mu ch-studied exponential function- 
als of Brownian motion (jYorl . Il992i or lYorL I2001L p. 15). In particular 
Dufresnel (1990) has shown that such a functional 



exp 



(« 



B) 



bs ] ds 



has the distribution of 2 / (a 2 T 2 b / a 2 ) i where T K is a Gamma-distributed 
random variable of index k. In summary, 

Theorem 3.1 Let T coup ii ng be the coupling time for two-dimensional 
Brownian motion plus Levy stochastic area under a mixture of reflec- 
tion coupling and synchronous coupling. Using the notation above, let 
min{l,a /W } be the proportion of reflection coupling. Scale T coup ii ng 
by the square of the ratio between initial areal difference Uq and initial 
spatial distance Vb; if 1 < a 2 < 2 then as Wo = Uo/V 2 — ?> 00 so a 
re-scaling of the coupling time has limiting Inverse Gamma distribution 



Y±\ T 

jj^ I -L coupling 



1 



2 + a 2 T 



(2-a 2 )/(2(2+a=)) 
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Thus the tail behaviour of T C oupiing is governed by the index of the 
Gamma random variable, which here cannot exceed i (the limiting case 
when a 2 — > 1). This index is unattainable by this means since H — 2K 



Kendahl (|2007l Section 3) exhibits a similar coupling for the planar 



case in which there is a state-dependent switching between reflection 
and synchronous coupling, depending on whether the ratio W exceeds a 
specified threshold; it would be interesting to calculate the tail-behaviour 
of the inverse of the coupling time in this case. 

Note that scaling by (Uo/Vo) 2 rather than Wq = Uo/Vq quanti- 
fies something which can be observed from detailed inspection of the 
stochastic differential system (|3.5p ; the rate of evolution of the areal dis- 
tance U is reduced if the spatial distance V is small. The requirement 
Wq — > cxd is present mainly to remove the effect of higher-order terms 
(dr/W 2 ) in systems such as p.9[) . and in particular to ensure in (|3.13p 
that min{W^ 2 ,a 2 } = a 2 for all time with high probability. 

In fact one can do markedly better than the reflection-synchronous 
mixture coupling of Theorem 13.11 by replacing synchronous coupling by 
a rotation coupling for which sin 9 = \f~2~ f3/W: similar calculations then 
show the index of the inverse of the limiting scaled coupling time can 
be increased up to the limit of ^, which remarkably is the index of the 
inverse of the coupling time for reflection coupling of Brownian motion 
alone! However this limit is not attainable by these mixture couplings, 
as it corresponds to a limiting case of (3 = 2, a 2 — 3. At this choice of 
parameter values H — 2K again behaves like a Brownian motion so we 
do not have W —> oo. 

In higher dimensions similar calculations can be carried out, but the 
geometry is more complicated; in the planar case the form of Z_ is essen- 
tially constant, whereas it will evolve stochastically in higher dimensions 
and relate non-trivially to v. This leads to correspondingly weaker res- 
ults: the inverse limiting scaled coupling time can be bounded above and 
below using Gamma distributions of different indices. 

We should not expect these mixture couplin gs to be maxim al, even 
within the class of co- adapted couplings. Indeed iKendalll (|2007l ) gives a 
heuristic argument to show that maximality amongst co-adapted coup- 
lings should be expected only when one Brownian differential is a (state- 
dependent) rotation or rotated reflection of the other. The interest of 
these mixture couplings lies in the ease with which one may derive lim- 
iting distributions for them, hence gaining a good perspective on how 
rapidly one may couple the stochastic area. 
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4 Conclusion 



After reviewing aspects of coupling theory, we have indicated an ap- 
proach to co-adapted coupling of Brownian motion and its stochastic 
areas, and shown how in the planar case one can use Dufresne's formula 
to derive asymptotics of coupling time distributions for suitable mixed 
couplings. Aspects of these asymptotic distributions indicate the price 
that is to be paid for coupling stochastic areas as well as the Brownian 
motions themselves; however it is clear that these mixed couplings should 
not be expected to be maximal amongst all co-adapted couplings. Ac- 
cordingly a very interesting direction for future research is to develop 
these methods to derive estimates for coupling time distributions for 
more efficie nt couplings usi ng state-dependent coupling strategies as ex- 
emplified in iKendalll (|2007l Section 3). Progress in this direction would 



deliv er probabilistic gradient estimat es in the manner o f lCranstonl (|1991 
(contrast the analytic work of iBakrv et al. 2008) 

A further challenge is to develop these techniques for higher- 
dimensional cases. Here the two-dimensional approach extends naively 
to deliver upper and lower bounding distributions; a more satisfactory 
answer with tighter bounds will require careful analysis of the evolution 
under the coupling of the geometry as expressed by the pair (i/,Z_). 

A major piece of unfinished business in this area is to determine the ex- 
tent to which these co-adapted coupling results extend to higher-order it- 
erated path integrals (simple Ito calculus demonstrates that it suffices to 
couple Levy stochastic areas in order to couple all possible non-iterated 
path integrals of the form J Bi dBj ) . Some tentative insight is offered by 
the role played by the Mo rse-Thue sequence for iterated time-integrals 
( Kendall and Price . 2004h . Moreover it is possible to generalize the in- 
variance considerations underlying (|3.4[) for the areal difference, so as 
to produce similarly invariant differences of hig her-orde r itera ted path 
integrals. But at present the closing question of IKendalll (|2007l ) still re- 
mains open, whether one can co-adaptively couple Brownian motions 
together with all possible iterated path and time-integrals up to a fixed 
order of iteration. 
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